Mathematical modeling of filamentous 
microorganisms 

Michele Bezzi* and Andrea Ciliberto^ 

* SONY CSL, 6, rue Amyot, 75005, Paris, France. 

t Dept. Biology, Virginia Polytechnic Institute, 
Blacksburg, VA, 24061, USA. 

Abstract 

Growth patterns generated by filamentous organisms (e.g. actino- 
mycetes and fungi) involve spatial and temporal dynamics at different 
length scales. Several mathematical models have been proposed in the 
last thirty years to address these specific dynamics. Phenomenologi- 
cal macroscopic models are able to reproduce the temporal dynamics 
of colony-related quantities (e.g. colony growth rate) but do not ex- 
plain the development of mycelial morphologies nor the single hyphal 
growth. Reaction-diffusion models are a bridge between macroscopic 
and microscopic worlds as they produce mean-field approximations of 
single-cell behaviors. Microscopic models describe intracellular events, 
such as branching, septation and translocation. Finally, completely 
discrete models, cellular automata, simulate the microscopic interac- 
tion among cells to reproduce emergent cooperative behaviors of large 
colonies. In this comment, we review a selection of models for each of 
these length scales, stressing their advantages and shortcomings. 

1 Introduction 

The panorama of mathematical models of filamentous microorganisms is 
quite large, for several reasons. To begin with, filamentous organisms grown 
on solid media, give rise to morphologies (circular colonies with branched 
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mycelia) that are easily observable but need to be put in the rigorous frame 
of mathematics to be properly understood (see [6, 5, 36] for examples of 
mathematical models of microorganisms growth patterns). Qualitative argu- 
ments that have proved very useful in other areas of biology simply could not 
match the dynamical and spatially extended problems introduced by mycelial 
growth. A second reason of interest concerns industrial applications of fil- 
amentous organisms, since they are among the main producers of enzymes 
and antibiotics (for a recent review on fungal morphology and industrial ap- 
plication see [42]). Not surprisingly, several mathematical models have been 
produced to improve growth and catabolyte production of fungi and acti- 
nomycetes. Finally, multicellular life-style of bacteria has recently attracted 
much attention, and particularly developmental processes related with mul- 
ticellularity in prokaryotes [53]. Among bacteria, Streptomocyes is one of the 
most striking examples of a multicellular bacteria, and differentiation is a 
well established and documented process in this genus. 

In this review, we will not cover extensively the published models of 
growth of filamentous organisms, rather we will discuss a few examples of 
models at different organization level, ranging from the single cell to the 
whole colony level. Modeling mycelial growth can be undertaken at different 
levels of complexity and length scale. Microscopic models are focused on 
single cell scale, so they deal with tip growth, branching and septation (see for 
example [14, 33]), while macroscopic models study the evolution of quantities 
on the scale of the whole colony, such as biomass or total hyphal length (see 
for example [56, 15, 16]). 1 

In principle, microscopic, single-cell description can be applied individ- 
ually to an huge number of cells, to simulate the macroscopic behavior of 
the whole colony. Unfortunately, this approach presents many shortcomings: 
first the many non-linearities present in these models require large scale sim- 
ulations which can be computationally inefficient and unfeasible even using 
modern computers. In addition, a big number of equations implies a big 
number of parameters, which are often difficult to estimate from experiments. 
The uncertainty of their values results in many different possible outcomes in 
simulations. Finally the results are often difficult to understand and analyze, 

1 Notice that "microscopic" and "macroscopic" in this context do not always correspond 
to the same terms used in other biological or physical modeling contexts, for instance in 
our case microscopic means the level of a single cell or a single hypha, which is in fact a 
macroscopic physical object itself, with its own complex dynamics. 
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so that even relatively simple phenomena can be obscured by the complexity 
of the model. 

A possible link between microscopic and macroscopic approaches con- 
sists in models whose variables are local densities of hyphae and concentra- 
tions of growth-related substances. In these models, which can be viewed 
as a local mean-field approximation of microscopic models, single-cell in- 
dividuality and local fluctuations in the number of cells or molecules are 
lost, nevertheless some spatial structure is retained so that they are able to 
reproduce some of the macroscopic growth patterns of the colony. Mathe- 
matically, these are reaction- diffusion models [19, 20, 21, 60, 9], composed 
of a set of coupled partial differential equations with reaction and diffusion 
terms. Reaction-diffusion models have been widely used, especially to de- 
scribe pattern-formation. 

In the following sections, we first will give a short summary of the main bi- 
ological features of growth of filamentous microorganisms (Section 2). Then, 
we will present examples of macroscopic, microscopic, reaction-diffusion and 
completely discrete (cellular automata) models for mycelial growth in fungal 
and bacterial colonies (Section 3), stressing the advantages and limitations 
of the different approaches. Finally, in Section 4, we will discuss pattern 
formation for aerial mycelium. General conclusions are drawn in the last 
section. 

2 Growth of filamentous microorganisms 

In bacteria that divide by binary cell division, mother and daughter cells 
separate at the end of the cell replication process, with a minimal spatial 
displacement: they basically pile up on each other. Cell colony increases in 
size, and exploits the nearest resources, until they have been depleted. The 
limits of this strategy are particularly evident for growth on solid (typically 
agar) media, where colonies are quite small and stop growth after a few 
days. Particularly, these cells are neither able to escape from a region of low 
nutrient concentration, nor able to avoid depletion of resources by restricting 
growth to a few specialized cells. In the contrary, a filamentous prokaryote 
such as Streptomyces is able to colonize a whole agar plate, even in minimal 
medium. How is that possible? 

In filamentous organisms, cells are organized in hyphae, filaments of cells 
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divided by septa (for reviews on filamentous organisms growth and life cycle 
see Refs. [48, 57]). At cell division, mother and daughter cells are not phys- 
ically separated, but are kept together by the septum, Fig. 1. This physical 
constraint has wide effects on the mechanisms of growth and division: are 
mother and daughter cells identical concerning growth and division? The 
answer to this question lies at the core of the strategies implemented by fila- 
mentous organisms to exploit environmental resources, and consequently to 
the growth kinetics and patterns observed on solid media. 

In this section, we set the stage to answer this key question, describing 
the basic mechanisms of cell division and hyphal growth, and we explore the 
most common morphologies and growth kinetics of filamentous organisms on 
solid media. 

2.1 Hyphal growth 

A typical colony of fungi or actinomycetes grown on solid media is formed 
by two different classes of mycelia: vegetative and aerial mycelium. The for- 
mer penetrates the medium, and is responsible for absorbing nutrients from 
the environment, while the latter, developed from the vegetative mycelium, 
releases spores in the air. 

We shall describe hyphal growth and kinetics following the development 
of a colony, from germination to the formation of a fully developed colony. 
A rigorous general description of filamentous organisms is impossible since 
growth mechanisms differ between fungi and actinomycetes, and even among 
one actinomycete (or fungus) species to another. However, it is possible to 
describe some general properties which apply to a hypothetical filamentous 
organism, keeping in mind that every specific case will somehow differ from 
this scheme. 

A colony originates from a single spore, which germinates in favorable 
conditions. Cell elongation is localized to the tips, where new cell wall is 
continuously assembled [44, 25]. In a typical cell division cycle, Fig. 1, 
the cell grows, replicates its DNA, and the genetic material is segregated to 
mother and daughter cells. Finally, once the mother cell has reached a critical 
size (roughly, it has doubled its original size), septum formation separates 
mother and daughter cells [30, 23]. 

This cell division process gives rise to an apical and a subapical cell. They 
differ greatly, since the apical cell is actively growing, growth being localized 
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again at the hyphal tip, while the subapical cell does not grow, but sup- 
ports tip growth by producing material that is delivered to the apical cell. 
Successive rounds of coordinated elongation and septa formation produce a 
hypha formed by several different cells separated by septa [30]. Notice the 
difference of this process from cell division occuring in cells that replicate by 
binary division: in a colony of E. coli every cell of the population can repli- 
cate. In filamentous organisms the very mechanism of cell division implies 
differentiation for cells located in different positions along an hypha. 

The hypothetical organism we have described so far is formed by a hypha 
growing at the tips. Extension rate increases exponentially during the first 
hours, synchronous with periodical rounds of DNA synthesis, and with the 
lengthening of the supporting hypha [31]. Tip growth rate reaches a maxi- 
mum when growth rate equals the rate of transport to the tips of material 
needed for cell wall construction. After this phase, tips growth cannot match 
the production of cell wall material and elongation rate is constant. How- 
ever, the excess material is not wasted, but used by a second growing point 
branched from the hypha, Fig. 1. The very same growth properties described 
for the original hypha apply to the new branch: cell growth is localized at 
the tips, new tips can be formed by branching. 

Generally, in rich media branching is highly favored, while it is inhibited 
in poor media. However, during early growth, both the total number of 
branches and the total mycelial length increase exponentially with the same 
specific growth rate (//), so that the colony expands exponentially. Their 
ratio is called the hyphal growth unit, and it reaches a constant value after 
an oscillatory transient ( a more quantitative analysis of this phase of growth 
is given in Section 3). 

The relationship between branching, growth and septa formation varies 
with the species and with the strains, both in actinomycetes and fungi. For 
example, Streptomyces granaticolor shows a pattern similar to fungi where 
branches are localized near to septa [30], while in Streptomyces coelicolor 
branching occurs far from them [48]. Moreover, in Streptomyces coelicolor 
apparently branching and septation can be completely separated, as mutants 
unable to septate are still capable of branching, and develop a vital mycelium 
[41]. 
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2.2 Colony morphology and kinetics 

After a few hours, a colony usually develops a circular shape, and the colony 
radius elongates with a linear extension rate, driven by the growing tips at 
the colony margins (see Fig. 2). More precisely, it is possible to define a 
"peripheral growth zone" with width w [55], actually a ring of hyphae, that 
is actively involved in colony elongation. Both hyphal tips, where the growth 
process occurs, and supporting cells, that produce material needed for growth 
to be delivered to the tips, belong to the peripheral growth zone. 

An empirical equation describes the relationship between colony radial 
rate (dr/dt = K r , at this stage, colony radius r is an obvious choice to 
measure colony growth), w and the specific growth rate \i [55] 

^ = K T = fj, r w (1) 

This empirical relationship is useful to understand the strategy adopted 
by filamentous organisms to exploit resources [48]. Both in a rich and poor 
environment, a colony grows to completely cover an agar dish, but the to- 
tal biomass produced is different in the two media. When the environment 
is poor in nutrients, branching is kept low, the number of tips is low and 
therefore w can be kept high. This way, even though \i is small, the colony 
moves with an high extension rate K r through a region that cannot pro- 
vide enough nutrients to maintain a large number of cells. That is why, for 
example, Streptomyces can fill a whole petri dish even growing in minimal 
media. On the other hand, if the area is rich in nutrients, branching is very 
high, to exploit completely the available resources. The resulting w is small, 
because many growing tips need to be supported, but the extension rate K T 
is kept high by a large fi. A very similar strategy is adopted by fungi, but 
the dimensions are different: w in Streptomyces are typically one order of 
magnitude smaller than in fungi. 

Even though only cells located at the margins are responsible for radial 
growth, growth occurs in the center of a colony as well. In Streptomyces gra- 
naticolor, Kretschmer [30] reports various cycles of DNA replication, branch- 
ing and septation in subapical cells, as certified by the smaller average size of 
subapical not branched cells compared to apical cells. However, eventually, 
branching and septation cease completely in the center of a colony. Cells 
located in this area are the first to detect a lack of nutrients and the pres- 
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ence of a high concentration of catabolites [52]. Not surprisingly some of 
them differentiate into aerial mycelium, that will give rise to spores. These 
specialized cells will be transported by physical or biological vectors to new 
regions where they will develop into new colonies. 

A detailed description of aerial mycelium formation is beyond the scope 
of this comment (see [12]); however, it is important to stress that differentia- 
tion relies upon intercellular signaling. For example, in S. griseus a substance 
termed A-factor is involved in triggering aerial mycelium production [26] . Ex- 
tracellular diffusion of the inducers of aerial mycelium guarantees a global 
coupling of the cell population. Aerial mycelium growth is supported by other 
cells that undergo lysis, thus providing nutrients for the aerial mycelium [40]. 
At this stage, Streptomyces produce antibiotics to prevent other bacteria to 
have access to this source of nutrients [11]. Usually aerial mycelium grows in 
a compact circular pattern but different patterns can develop as well. S. rut- 
gersensis grown on minimal media forms concentric rings of aerial mycelium 
(see Fig. 7) [7]. Similar but order of magnitude larger rings are also ob- 
served in Neurospora crassa [17]. Mathematical models have been produced 
to explain the generation of these patterns: they will be presented in Section 
4. 

3 Modeling fungal growth 
3.1 Macroscopic models 

One of the first empirical macroscopic model was proposed by Trinci [56] to 
explain mycelial growth in its early stage. Let us call ip(t) the number of 
branches and X(t) the total mycelial length at time t. During the early stages 
of growth, the temporal evolution of these quantities is well reproduced by 
the following equations: 



m 
\(t) 



aexp(/it) 
cexp(/it) 



(2) 
(3) 
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with a and c constants and /i the growth rate. From Eqs. (2-3), we immedi- 
ately get: 

<«<) _ r 

that is, the ratio G (called hyphal growth unit) is constant over time. 

These phenomenological relations are in agreement with growth patterns 
of several filamentous fungi (see for example [10]) and bacteria (see for ex- 
ample, Figure 2 in [2]). 

During the initial phase of growth, after a short lag period, also the radius 
of the colony r c grows exponentially, (see Ref. [29] for a model of this phase of 
growth). Then, as described in Section 2, linear radial growth sets up while 
hyphae within the interior of the colony do not grow further [56]. During 
this phase, there is a quadratic increase in biomass. In fact, 

r c oc t 

(with t the time from the beginning of linear growth phase) and assum- 
ing biomass M uniformly distributed throughout the colony, (i.e. the total 
biomass proportional to the colony surface). 

M oc r 2 M (xt 2 

That is, total biomass grows quadratically in time. 

This simple model found many experimental confirmations, but it does 
not supply any information on the local distribution of the biomass, and 
in particular on the different behaviors between center and border of the 
colony. Using a reaction-diffusion model Ferret et al. [22] have predicted the 
evolution of the biomass and tips density under different growing conditions. 
In presence of unlimited nutrients and no spatial constraints, Eqs. (3-2) are 
recovered, and the biomass increases exponentially with time. However, as 
mentioned above, growth stops in the center of the colony: in the model the 
same happens when the biomass density m(x, t) reaches a maximum value, 
n^max, while at the colony border the growth continues freely. These differ- 
ential growing velocities are modeled by introducing a collision probability P c 
as a function of local biomass density m(x,t). In the the central part of the 
colony, the local density quickly reaches its maximum value, m(x, t) = m max , 
there is no free space available for growth, collisions are unavoidable and the 
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free space (proportional to 1—P C ) is zero: the mycelium cannot grow further. 
On the other hand, at the colony border, no collisions happen, P c = 0, and 
the hyphae can grow freely outward. Assuming a linear model for P(m(x, t)) 
(see Fig. 2 of the cited paper), that fulfills these two limiting cases, the model 
is able to predict the temporal behavior of total biomass and the colony bor- 
der evolution (i.e. internal and external radius, defined as the radius of the 
not-growing core, and the radius of the whole colony, respectively). The 
model has been validated by comparing its predictions to the behaviors of 
two different fungi: Gibberella fujikuroi and Aspergillus oryzae. 

These macroscopic models are suitable for modeling the growth of the 
whole colony, but they are not able to grab all the complex fine structure of 
mycelial growth. In fact, in addiction to tip extension and branching, filamen- 
tous growth is characterized by septation and internal transport mechanisms 
(translocation): these microscopic features will be taken into account in the 
next class of models. 

3.2 Microscopic models 

Septation and branching Once the nuclear material is replicated, mycelial 
microorganisms divide a hyphal compartment in two parts by septation, 
Fig. 1. Branching usually occurs after septation, with a fixed time delay 
(see for example the case of Geotrichum candidum in [23]). A model that 
includes branching and septation has been introduced by Yang et al. [60] The 
model includes a deterministic part, that mirrors tip growth and septation, 
and a stochastic part to keep track of new branch initiations. Septation is 
triggered when the amount of nuclear material has been doubled and segre- 
gated. Due to the random behavior of branching, new branches are supposed 
to be generated at random around their respective septa according to a nor- 
mally distributed probability with a fixed variance a\ b . In the limiting case, 
o~\ = 0, new branches grow just behind septa. Tip growth and new branch- 
ing directions are calculated according to a random model described in [61]. 
In short, experimental observations in 2-D colonies reveal that tip growth 
angles Aj g and branching outgrowth angles Ajb are normally distributed 
around zero. In the model, that represents growth in a 3-D environment. 
The orientation of the primary hyphal segment with respect to the straight 
direction is obtained sampling the experimental A'-fg distribution. In the 
plane perpendicular to the growth axis, the direction is fixed by an addi- 
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tional random angle assumed to be uniformly distributed between [0, 2ir]. 
A simile procedure is used for determining the branching angle A7& in the 
plane perpendicular to the hyphal growth axis (see Fig. 4). 

Simulation results are in good agreement with experimental data for the 
fungus Geotrichum candidum (taken from Ref. [23]) and for the filamentous 
bacterium Streptomyces tendae. Both germ tube growth and temporal be- 
havior of the number of tips and septa can be reproduced. In addition, the 
morphology of the growing mycelium appears similar to that observed in 
experiments (see Fig. 10 in [23]). 

Translocation Another peculiar aspect of the internal kinetics of a grow- 
ing mycelium is the capability to distribute nutrients through the mycelial 
network (translocation), both passively (due to diffusion through the hyphae) 
and actively (due to active and hence energy-dependent metabolism) [45, 46]. 
This way, fungi can absorbe nutrients in some parts of the mycelium and 
transport them through the network, particularly where food is scarce. This 
mechanism can play a crucial role for colonies growing in highly heteroge- 
neous environments, such as most of the natural ones. 

Among those that have been recently proposed [15, 16], we will describe 
here a recent model introduced by Boswell et al. [9] to study the fungus 
Rhizocotonia soloni. The model includes five components (densities): active 
and inactive hyphal densities, tip density, internal and external substrate 
concentrations. Inactive mycelia, usually the center of the colony, do not 
contribute to translocation, while translocation occurs in active hyphae. In- 
ternal passive translocation is modeled as a diffusion- like term for internal 
substrate, s, that depends on mycelial concentration m a , i.e. in presence of a 
dense mycelial network the diffusion process is much faster than in a sparse 
one. Active translocation is modeled by assuming that active mechanisms 
are due to the continuous substrate demand by growing tips, i.e. the flux fol- 
lows the gradient of tips density pt- So, active translocation moves substrate 
from low-density tip areas to high density ones (so called tip-driven diffusion) 
through the mycelial network (i.e. it also depends on m a ). In mathematical 
terms: 



The parameters are obtained from experiments. Using this model the 
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authors explore the interplay between active and passive translocation in 
a growing colony in two kinds of growth media, homogeneous (nutrients 
are homogeneously spread on the dish) 2 and heterogeneous (substrates are 
patchy distributed). Simulations suggest that diffusive translocation is used 
for random exploration while tip-driven translocation is used for optimizing 
the distribution of available resources. This view is coherent with the idea 
that optimal search combines a random process with some more complex, 
but usually more demanding, strategies of resource exploitation. 

Single hypha models Going down in the length scale, we find single 
hypha models. The average kinetics implemented in the previous models 
do not describe properly the evolution of individual hyphae. Indeed, due 
to the complexity of the hyphal structure and the different growth rates of 
individuals branches within an hypha (see [56]), stochasticity has to be taken 
into account. 

Christiansen et al. [14], studied the growth of a single hyphal element of 
Aspergillus oryzae, combining image analysis with a model that makes use 
of Monte Carlo simulations. They consider the growth starting from a single 
spore; at each time step (5r ), the probability of branching is: exp(gf, ran 5r) — 
1. The branching rate qv an is computed from a phenomeno logical model 
of tip kinetics qu P 3 , and from observed distribution of the number of new 
branches as function of hyphal length n\. From experimental observations, 
tip velocity as function of hyphal length I can be fitted by the following 
function (see Fig. 3): 

Qu P (l) = ktipj—^r (5) 

where ku p is the maximum tip extension rate, K s is an empirical saturation 
constant, that depends on environmental conditions (glucose concentrations), 
and varies among different types of branches (K s is usually larger for the 
germtube than for secondary hyphae). 
In the model K s is calculated from: 

K s = exp((3p + a) 

2 Thc effects of different, but homogeneous, growth media have also been studied in 
Aspergillus Oryzae colonies [35, 38, 39], and will be discussed further in the case of Strep- 
tomyces colonies in Section 4. 

3 A detailed model of the elongation process of the growing tip can be found in Rcf. [51]. 
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with p a normally distributed random number, and the parameters a and (5 
fitted from experimental values of K s . Eq. (5) is a good approximation of the 
average behavior of tip growth rate q tip ] but detailed experimental measure- 
ments show that a decrease in tip extension rate occurs in correspondence 
to the formation of new apical branches (see Fig. 3 adapted from Fig. 2. in 
Ref. [14], but also the discussion in Ref. [28] where no correlations between 
tip extension rates and branching are found). Following Ref. [14], we will 
neglect these fluctuations in the current analysis. After the hyphal length 
has reached a certain threshold l m i n (that depends on the individual hypha, 
and differs greatly between primary, germtube, and secondary hyphae), the 
number of new branches n\ formed on a hyphae grows linearly as function of 
length: 

{0 . I < lmin / r \ 
Ul - I ■ ) ■ I > I ■ (6) 
"^V 1 L min) ■ 1 ^ l rmn 

where kb is a branching constant (actually, it assumes two different values 
for the germtube and for secondary hyphae). From Eqs. (5) and (6), the 
branching intensity (i.e. the number of branches per time unity) can be 
derived: 

{0 . / < l m in /_x 
Ks+l ' l min 

Using these values of qbran we can simulate the growth of a single hyphal 
element using a standard Monte Carlo technique. In a standard simulation, 
a single spore has been considered, the tip grows with a rate q t i P and at each 
time step St a new branch is formed with a probability exp{qbranSr) — 1. 
Although hyphae grow only at the tips, with a saturating velocity, the total 
length of mycelium increases exponentially due to the exponential increase 
in the number of branches (after the minimal length l min has been reached). 
The simulated growth kinetics agree with the observed experimental patterns 
for the total length and the number of tips for a single hypha. 

A detailed study of single hypha growth has been presented as well by 
Hutchinson et al. [28] . They record the temporal evolution of primary hyphae 
and first and second order branches in six colonies a few hours after germina- 
tion. Using the experimentally observed distributions, they proposed a model 
of mycelial growth: hyphae perform a straight-line growth for a length sam- 
pled from an experimental fitted distribution. At this point, branching can 
occur, and the direction for the new branch is obtained by sampling from 



12 



experimental distributions. This model is capable of reproducing the mor- 
phology of the early stages of Mucor hiemalis growth, and represents one 
of the first model where a single hypha dynamics (microscopic) is used for 
generating macroscopic growth patterns (morphology of the mycelia). 

From the single hypha level, we could go down one level in the analysis 
to the genetic and molecular basis of filamentous growth (i.e. to the genetic 
networks underlying these processes). However, although progresses have 
been recently made in elucidating the molecular basis of mycelial growth 
and development (for example, see [13] for a review on Streptomyces) , no 
mathematical link between physiology of cell behavior and molecular inter- 
actions have been successfully produced so far. 

3.3 Cellular Automata 

Cellular Automata (CA) are fully discrete dynamical systems, i.e. their dy- 
namical variables are defined at the nodes of a lattice (spatial discretization), 
with values taken from a finite set (state discretization) and evolve in discrete 
time steps (temporal discretization). They were introduced more than half 
century ago by John Von Neumann's [59] . As the choice of the name suggests, 
the original idea is that the complex behavior characterizing many biological 
systems should originate as a collective effect from simple local interaction 
rules (cell- interactions). Since then, CA have been applied to a vast range of 
biological and physical problems (see for example the reviews [18, 8, 3] and 
the "classical" works [59, 24, 62]) 4 

A CA model for fungal growth has been proposed by Lazlo & Silman [33]. 
They consider a two dimensional square grid (that mirrors the dish) where 
each site has 8 neighbor sites (Moore neighborhood). Each site can be either 
occupied by a mycelium or empty. A probabilistic rule is set for updating the 
status of each site of the lattice simultaneously. In its simplest version the 

4 Another class of fully discrete models that emphasize local interactions as opposed 
to global ones are Lindenmayer Systems (L-systems) [34, 50]. L-systems are basically 
a set of string re-writing rules applied recursively. They are closely related to cellular 
automata, but usually they lack an explicit spatial extension. These rules can mirror 
growth process, and they are generally used for generation of fractals and modeling of 
plants. An applications to fungal growth, Mucor hiemalis, was proposed by Soddcl et 
al. [54]. We will not discuss these models here, but we refer the readers to the cited paper 
and references therein. 
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rule is the following: each empty site will become occupied with a probability 
p a dependent on the number of neighbor occupied sites n£. 

• rid = 1 (one neighbor site occupied) =>- p a = 0.12 

• n d = 2 Pa = 0.25 

• n d = 3 => p a = 0.50 

• ™d > 4 (overcrowded) => p a = 

Even using only this simple rule, interesting patters are generated, and 
the temporal behaviors of biomass and radius of the colony are well repro- 
duced. This basic model has been extended considering three different states 
of differentiation: vegetative mycelium, differentiating mycelia (an interme- 
diary state) and spore formation, and setting the corresponding probability 
of transition according to local occupancy rules. Still keeping the model 
simple (one single byte is enough to represent the state of a site) the au- 
tomaton is able to produce patterns for spore and mycelia, reminiscent of 
the experimental growth patterns. 

One of the advantages of CA, is that they can be used in testing hy- 
potheses about the interactions among cells: some schematization (rules) of 
microscopic dynamics can be easily implemented as cellular automata; then 
we can observe the macroscopic patterns generated and compare it with ex- 
perimental patterns. This class of models is particularly useful when the 
macroscopic patterns observed do not depend on the fine details of micro- 
scopic interactions, but they emerge as a collective behavior of many cells 
(the site of the automaton). In this framework, cellular automata, despite 
their simplicity, turn out to be a valuable tools for the studies of emergent 
cooperative effects. 

4 Aerial mycelium pattern formation 

As described in Section 2.2, aerial mycelium can form concentric ring patterns 
in S. rutgersensis colonies (see Fig. 7). Ring patterns are found also in other 
bacterial and fungi colonies [1]. Many models can produce such patterns [37]; 
in the case of bacterial colonies, it was suggested that the interplay of front 
propagation and Turing instability leads to concentric rings and spot patterns 
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[58]. In the more specific case of filamentous bacteria, a different approach 
based on competition for resources has been proposed by Bezzi et al. [7]. The 
model is based on the idea that growing hyphae have the biological objective 
to find nutrients to give rise to spores (aerial mycelia). As a result, we expect 
a strong competition to arise on minimal media for the energetic resources 
between neighbor substrate mycelia, whereas in maximal media, where there 
are sufficient nutrients, the competition is weaker. 

If cells are connected mainly along the radial direction (hyphae start 
growing mainly radially from the initial spore), then competition will be 
stronger along this direction than along the tangential one. In other words, 
in the growing edge of the colony, the competition is not isotropic but, fol- 
lowing the vegetative mycelium morphology, it will be stronger among cells 
belonging to neighboring circumferences (radial direction) than among cells 
belonging to the same (tangential direction). 

In presence of high concentrations of nutrients, the competition is rather 
weak and the final distribution of aerial mycelium has the shape of a gaussian- 
like distribution centered on the initial spore (similar to the one shown in 
Fig. 5, right). If the colony grows in a nutrient-limited environment, hyphae 
that are neighbor along the radial direction compete strongly for nutrients, 
aerial mycelia development is inhibited, and concentric ring patterns appear 
(see Fig. 5, left). Despite its simplicity (one differential equation) the model is 
able to capture many peculiar features of Streptomyces rutgersensis growth. 
However, since the model does not include explicitly the growth of vegetative 
mycelium, some important aspects, like the "interference" patterns among 
neighboring colonies (Fig. 6), are not reproduced. Clearly, a more detailed 
view is necessary; a coupled maps model including vegetative mycelia and 
tip growth is under development [4] whose results are shortly summarized in 
Fig. 8. 

As in the case of vegetative mycelia growth, different approaches are re- 
quired to describe the microscopic details, in which fluctuations and stochas- 
ticity play a major role. Monte Carlo simulations are a possible choice, an 
alternative strategy is using cellular automata, in which simple microscopic 
local rules (deterministic or probabilistic) can produce complex patterns. 

Cellular automata have also been used for studying pattern formation 
in bacteria, in particular for Streptomyces rutgersensis [43]. In this two- 
dimensional model the building block is the individual cell. A hypha is com- 
posed of several cells, which do not move in space. Cells absorb food (energy) 
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according to its local concentration, convert and store it internally and use 
the stored energy for metabolism, growth and reproduction. If the amount 
of stored energy lowers below a given threshold, the cells start to starve. The 
model assumes that this starvation phase triggers the formation of aerial 
mycelium. When a cell reaches a given dimension and has accumulated suf- 
ficient energy, it reproduces by binary division or gemmation. Despite many 
approximations introduced, the resulting patterns of vegetative mycelium are 
very reminiscent of experimental patterns: at the beginning of the life of the 
colony (i.e., for the first generations) the hyphae are distributed randomly in 
space, but then, when food becomes scarce, they are selected to align with 
food gradient. 

Similar concentric rings patterns are also observed in fungi, among them: 
Neurospora crassa. In Ref. [17], Deutsch et al. present a cellular automaton 
for modeling Neurospora crassa growth patterns. The two-dimensional pat- 
terns are reproduced using a one-dimensional cellular automaton; but at each 
time step the one dimensional grid (called "base space" for the automaton) 
represents a circle of increasing radius (so that the radial structure can be 
preserved). The model is capable of reproducing the morphologies of a single 
growing colony of Neurospora Crassa, but, due to the peculiar implementa- 
tion of the radial symmetry, it is not applicable to the studies of interference 
patterns (similar to S. rutgersensis pattern shown in Fig. 6). 

5 Conclusions 

In filamentous organisms (e.g. fungi and actinomycetes), mother and daugh- 
ter cells never detach: although they are separated by a septum, they remain 
physically connected. A filament of cells separated by septa is called an hy- 
pha, which is the building block of filamentous organisms. Hyphae are usu- 
ally branched, and the complex of hyphae forming a fungal or actynomicete 
colony is called mycelium. 

The mycelial organization has deep consequences on growth of filamen- 
tous organisms, because active growth is localized only at hyphal tips. Colony 
growth rate is not simply dictated by hyphal tips growth, but also by the 
frequency of formation of new tips - branching frequency. Spatial organiza- 
tion and temporal growth are intrinsically connected, and the morphology of 
a colony is a key element to be taken into account if we want to understand 
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the growth of filamentous organisms. 

Growth and mycelial pattern formation involve events occurring at very 
different length scales, ranging from a colony scale (macroscopic) to a single- 
cell scale (microscopic). For example, experimentally spatial organization 
is described by quantities like the angle between an hypha and a related 
branch, the distance between two consecutive branches, the distance between 
two consecutive septa, the time course of the total amount of branches, the 
frequency of branches along an hypha, and so on. Evidently, these quantities 
are measured at the microscopic scale, but the resulting mycelial morphology 
is observed at the colony scale. On the other hand, the temporal dynamics of 
colony growth (e.g. colony radius growth rate) is a macroscopic phenomenon, 
which is originated by tips growth, a single-cell event. 

Mathematical models are needed to describe such complex phenomena 
that involve spatial and temporal dynamics, microscopic and macroscopic 
experimental observations. Indeed, several models have been developed in 
the last thirty years to reproduce growth kinetics and mycelial patterns in fil- 
amentous organisms. Different biological questions generate different type of 
modeling, that concur to draw a complete picture of filamentous organisms 
growth and pattern formation. Such models differ for the different organ- 
isms they apply to, as well as for the different scales of the experimental 
observations they simulate. 

In this comment we review a selection of models, organized according to 
the scale of the experimental observation from which the model originates. 
We start with macroscopic models that reproduce colony growth rate experi- 
ments, then, moving down in length scale, to some reaction diffusion models, 
that include more detailed phenomena like septation and translocation, and 
finally we terminate with single hypha models, that concern single hypha 
growth and branching. We also review cellular automata models that imple- 
ment microscopic interaction but at the same time they are able to reveal 
macroscopic patterns. 

The last step of the length scale should be the molecular level, but, to 
our knowledge, there are no models that use the genetic data produced in 
the last twenty years to address growth kinetics and pattern formation in 
filamentous organisms. Although genetic analysis revealed some of the genes 
involved differentiation and the way they are connected (i.e. the underlying 
gene networks), models based on these data are still missing. Such models 
have been produced for several basic processes like the cell cycle, circadian 
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rhythm and signal transduction pathways (for a review see [27]). By using 
similar mathematical analysis it should be possible to develop similar models 
to understand and modify in a predictable way filamentous organisms. 
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Figure 1: Filamentous organisms are organized in hyphae, where single cells 
are separated by septa and growth is localized at the hyphal tips. As apical 
cells grow, DNA is replicated, and eventually a septum separates mother and 
daughter cells. Occasionally, subapical cells can undergo a branching process. 
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Figure 2: Vegetative mycelium at the colony border of S. rutgersensis: notice 
hyphal tips growing outwards. 
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Figure 3: Tip extension rate for a single hypha, see Eq. (5). Experimental 
points are spread around the solid line, with temporary decreases corrispond- 
ing to branching point (adapted from Fig. 2 in Ref. [14]). 
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Figure 4: Schematic representation of growth direction angle and branching 
direction angle. Adaptation from [61]. 
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Figure 5: Simulations of aerial pattern formation of S. rutgersensis. Left, 
growth on minimal media, right, growth on maximal media. From [7] 




Figure 6: Colonies of S. rutgersensis develop interference pattern when grown 
in proximity 
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Figure 7: S. rutgersensis growth in minimal (top) and maximal culture me 
dia (bottom, 3 colonies are clearly visible). Aerial micelium: white spots 
vegetative micelium not visible. 
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Figure 8: (Top) Temporal evolution in minimal culture media Bottom: Tem- 
poral evolution in maximal culture media The model is based on a coupled- 
map lattice (a dynamical system discrete in time and space), with four "pop- 
ulations": concentration of vegetative mycelium and tips, concentration of 
food in the dish, concentration of food absorbed. These populations evolve 
in two different timescales: a fast one, that includes tip growth, diffusion 
and translocation; and a slower timescale for micelium duplication. The 
combined signal of food scarcity and high absorbed food triggers stochastic 
aerial micelium production (plotted in the figures). 
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